A divergent correlation length in off-equilibrium glasses 
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In off-equilibrium dynamics we define a dynamical correlation length which is proportional to the 
size of the region in which the atoms move in a correlated way. General arguments indicate that 
this dynamical correlation length diverges at large times in the glassy phase. Numerical simulations 
for binary mixtures point toward the correctness of this prediction. 



The glass transition has many puzzling features, one 
being the presence of a divergent time scale without a cor- 
responding divergence of the correlation length as mea- 
sured by considering equilibrium correlations. In this let- 
ter we will show that a divergent correlation length may 
be identified in the glassy phase, below the glass tran- 
sition point, if we study the off-equilibrium correlation 
functions. 

The basic idea is quite simple. In a glass atoms are 
usually frozen in some position. If they try to move, 
the other atoms push them in the original position. This 
cage effect is at the origine of the glass transition. At low 
temperatures, but still in the liquid phase, the viscosity 
is finite (albeit very large) and some movements are still 
allowed. In many theoretical approaches jjj these move- 
ments correspond to the formation of relatively large flu- 
idified domains. All the particles in the domain move of 
a certain amount and the radius of the domain diverges 
when we approach the transition. The Vogel-Fulcher law 
may be obtained if the volume of these fluidified domains 
diverges at (T — Tk)^ 1 and the process is controlled by 
the energy activation which should be proportional to 
the volume. The existence of processes involving a large 
number of particles in the liquid phase near the glass 
transition has been recently shown in Q . 

When we decrease the temperature, the number of 
particles involved increases; these processes involve the 
crossing of higher and higher barriers and they are ex- 
tremely hard to observe in equilibrium simulations 
Here we study similar processes which are present when 
we quench the system from an high temperature to a tem- 
perature well below T g Q-q|. The dynamics of systems 
going toward equilibrium has been the subject of a wide 
theoretical interests 0. In the high temperature phase 
the energy usually approaches equilibrium exponentially 
fast. On the contrary in the low temperature phase the 
energy has corrections which decrease as a power of of 
time: 

E(t) = E x + A t-^ T l (1) 

The dependence of A(T) on the temperature is quite in- 
structive. There are two simple possibilities. In the first 
one, the approach to equilibrium is dominated by en- 
tropic barriers and the value of the exponent A(T) is 
roughly independent from the the temperature. In the 



second case the approach to equilibrium is dominated by 
energy barriers: the value of A(T) strongly depends on 
the temperature and vanishes linearly with the temper- 
ature. Depending on the nature of the system, the ap- 
proach to equilibrium may be dominated by entropic or 
by energetic barriers. Two well studied cases, which cor- 
respond to the two different behaviours, are the spinodal 
decomposition (|] and the spin glasses |J. 

If coming from the high temperature phase we quench a 
ferromagnetic Ising system at temperatures much smaller 
than the critical one, we find bubbles in which the spins 
have a similar value and the radius of these bubbles in- 
creases at i 1 / 2 with z = 2 ||. The radius of these bubbles 
coincide with our dynamical correlation length. The per- 
sistence time of a bubble of size R is proportional to R z : 
a bubble which is present at time t after the quench takes 
a time O(i) to disappear. Here A(T) = 1/2 = l/z. All 
exponents are temperature independent. The spin-spin 
correlation C satisfies the scaling form 

C(t,t Wl r) ee (a(x,t w )a{y,t)) ~S{^-,-^—) (2) 

where \x — y\ = r and R(t w ) oc tlJ z . Eq. (|^) implies aging 
JlfH : at fixed r, at t and t w both large, the correlation 
function depends only on the ratio t/t w . 

An opposite behaviour is seen in spin glasses. The lo- 
cal magnetizations (m(i) = i being a point of the 
lattice) fluctuate from point to point with zero average 
(at zero external magnetic field). In order to observe 
an interesting behaviour it is convenient to consider two 
copies of the same system (<r and r) and to introduce 
the overlap q(i) = a(i)r(i). The quantity q plays the 
same role of the magnetization in ferromagnets: the sus- 
ceptibility associated to q diverges at the critical point; 
the statistical expectation value of qi(i) is equal to m(i) 2 
and it may different from zero only below the transition. 
The dynamical correlation function of q gives important 
information on the dynamics. In this case the theoretical 
understanding is not as good as in the case of ferromag- 
nets; however intensive numerical simulations (^] show 
that (in three dimensions) the dynamic correlation length 
diverges as a power of the time and that the ^-correlation 
function has a behaviour which is qualitatively similar to 
that of a ferromagnet: 
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C q (t, t w ,x) = (q(x, t w )q(y, t)) « r~ u S q (— , ), (3) 

where also here R(t w ) oc ti, . The main differences 
among spin glasses and ferromagnets are: (a) the pres- 
ence of the prefactor r~ u (u> w .5); (b) the value of the 
exponent z, which is approximately given z c T/T c (z c w 7 
and T c is the critical temperature). Also the relation 
among z(T) and A(T) is different from that of ferromag- 
nets: the energy approaches the equilibrium value with 
an exponent X(T) ~ 2.5z(T) , which it is proportional 
to the temperature. 

Summarizing ferromagnets and spin glasses show a 
clear cut behaviour characterized respectively by entropic 
and energetic barriers. Which are the theoretical expec- 
tations for structural glasses? We can answer to this 
question if we assume that the mean field for some gen- 
eralized spin glasses can be applied also to structural 
glasses pr| , p^[ . This conjecture can be considered as a 
rationalisation of the Gibbs Di Marzio approach; as a 
byproduct it implies the correctness of the mode cou- 
pling theory in an appropriate time-temperature window 
& 

In this approach the decay of the energy has a rather 
interesting behaviour. In the mean field theory there 
are two time windows. In the first one, at relative short 
times, the approach to equilibrium is dominated by en- 
tropic barriers and the value of the exponent A(T) is 
roughly independent from the the temperature. At larger 
times the approach to equilibrium is dominated by cross- 
ing of energy barriers; in this regime the time evolution 
strongly depends on the temperature. Numerical simu- 
lations (done both with Hamiltonian and dissipative dy- 
namics [|]|]]) for structural glasses are consistent with 
this picture. Here we will show that in the first (entropy 
dominated) regime there is dynamic correlation length 
which increases as power of time. 

We present the results of a simulation for binary flu- 
ids. We consider a mixture of soft particles of different 
sizes. Half of the particles are of type A, half of type B 
and the interaction among the particles is given by the 
Hamiltonian: 

ff = £((d( ? ;) + d(fc)) 12 |x l -x fc r 12 , (4) 

i<k 

where the radius (d) depends on the type of parti- 
cles. This model has been carefully studied in the past 
jl4],[|jl|] . The choice ds/d,A = 1.2 strongly inhibits crys- 
tallisation and the system goes into a glassy phase when 
it is cooled. Using the same conventions of previous inves- 
tigators we consider particles of average radius 1 at unit 
density. It is usual to introduce the quantity T = /3 4 . For 
quenching from T — oo the glass transition is known to 
happen around T c = 1.45 iQ. For computational rea- 
sons we have slightly modified this model by introducing 
a cutoff, i.e. we have put to zero the interaction when 



Xk\ 2 > 3. The differences in various thermodynamic 
quantities with the original model are of the order of 1%. 
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FIG. 1. The function r(g(r) — 1) versus r, the values of the 
fit parameters (according to eq. (5)) are for a = 0: r g = 1.29, 
d = .89, (f> = —1.15. Superimposed there is an indistinguish- 
able fit with a = 1. 

Our simulation are done using a Monte Carlo algo- 
rithm, which is more easy to deal with than molecular 
dynamics, if we change the temperature in an abrupt way. 
Each particle is shifted by a random amount at each step, 
and the size of the shift is fixed by the condition that the 
average acceptance rate is about .4. Particles are placed 
in a cubic box with periodic boundary conditions and at 
the end of each Monte Carlo sweep all the particles are 
shifted of the same vector in order to keep the center of 
mass fixed (3). We have done simulations for system of 
many sizes. Here we present the results for systems with 
27000 particles, which correspond to a box of size 30. We 
need to use large systems in order to avoid finite volume 
effects (the size of the system must be much larger 
than the dynamical correlation length we study). The 
value of r is 1.8, which is deep in the glassy region (it 
corresponds to a temperature three times smaller that 
the transition temperature 

Let us recall the results for the density-density corre- 
lation function (g(r)) at equal time at equilibrium after 
a rapid quench. The correlation function is shown in 
fig.(|l]) at distances greater than 1.6; in this region it can 
be fitted as 

g{r) = 1 + Ar~ a exp(-r/R g ) sin{2nr/d + 0), (5) 

which correspond to a pair of complex singularities in 
momentum space (simple poles for a — I). The fitted 
value of R g strongly depends on a; we find r g — 1.29 and 
2.13 respectively for a = and a = 1. A fit done only 
at large distance (in the region r > 4) is better and gives 
similar values of R g (1.88 and 1.41). The value of the 
correlation length often strongly depends on the power 
in the prefactor and, unless one has extremely good data, 
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it is difficult to determine both the power of the prefactor 
and the rate of the exponential decay. 
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FIG. 2. The overlap q as function of s = t/t w for 
t w = 2 9 , 2 13 and 2 17 , averaged over 3 initial conditions, the 
errors being smaller than the points. 

In order to verify the phenomenon of aging we have we 
have introduced the quantity q(t w , t) [BJ defined as 

q(t w ,t) = ^ 

i 

Qi(t w ,t) = y^w(xi(t + t w ) -x k (t w )), (6) 
k 

where the sum over k is done over particles of the same 
type of i. We have chosen the function w in such a way 
that it is very small when x » a and it is near to I 



for x < a, i.e. w(x) 
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i 12 ), with a = .22. 



The value of q will thus be a number very near to 1 for 
similar configurations (in which the particles have moved 
of less than a) and it will be much smaller (O(10 -1 )) for 
unrelated configurations; using the terminology of spin 
glasses q can be called the overlap of the two configura- 
tions. In fig. ^ we plot the overlap as function of s = t/t w 
at different values of the waiting time, i.e.t w — 2 9 , 2 13 
and 2 17 . The data weakly depend on the waiting time. It 
is possible to assume that q goes to limit when t w — » oo 
at fixed s, this limit being reached from above at small 
e and from below at large e. This way of approaching 
the asymptotic limit is quite a common feature in other 
systems, e.g. spin glasses fl7| . 

Our aim is to study the correlations of the particles 
which have moved in a sizable way when the time changes 
form t w to s t w . These particles, unless they have ex- 
actly replaced other particles, contribute to the decrease 
of the value of q and have a correlation in in position 
space which strongly depends on waiting time. To ev- 
idenziate this effect H||, for each particle and pairs of 
configuration at time t w and s t w we define the quantity: 



In this way J^i a i = 0- When the average value of q is 
not far 1/2, this procedure corresponds to put <Xj w 1, 
if the particle has moved less than a, and <Xj as —1, if 
the particle has moved more than a. Movements which 
correspond to an interchange of particles of the same type 
have no effect on a. 
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FIG. 3. The correlation c(r,t w ) as function of the distance 
for t„ = 2 9 , 2 13 and 2 17 

We now we define the function f(r,t w ) as the correla- 
tion of the particles at time s t w , where the contribution 
of two particle is weighted by a factor OiCJk- Equivalcntly, 
we can define the function = ^ i 5(x — Xi(t w ))(Ji(t w ), 
which is the local density p(x) multiplied by a. In this 
way we can write f(r) — (fi(x) [i(y)) , where \x — y\ = r. 

The correlation f(r) go to zero at large distances, by 
construction. We expect also that the values of r at which 
fir) is sizable different from zero correspond to the values 
of distance of particles which move in a correlated way. 

The results for the correlation f(r,t w ) computed with 
s — 3 are shown in fig. ^, where we plot c(r,t w ) = 
f(r,t w )/g(r) (we have divided the data for f(r,t w ) by 
g(r) in order to eliminate a natural oscillatory effect). 
At short times the correlations are present only at short 
distances, when the time increases they extend on a much 
larger region. We have analyzed the data by fitting 
c(r,t w ) as 



c(r,t w ) = exp(-r/R(t w ))(ci 



C2T 1 sin(27rr / d + <f)) , 
(8) 



(7) 



where we have taken the value of d from the previous fit 
of the correlation function g(r). (The same qualitative 
dependence of the dynamic correlation distance R(t w ) is 
obtained also using other fitting procedures; the precise 
value of R(t w ) being also here quite sensitive to the in- 
troduction of a power decaying prefactor). 

In fig. U we show the dynamic correlation distance as 
function of time. We see that the behaviour of 7Z(t w ) 
is well approximated by a power of time, the exponent 
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being about 0.16. 
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FIG. 4. On a double logarithmic scale we find the values 
of the correlation functions (R(t m ), axis on the left) as func- 
tion of the waiting time and the values of the integral of the 
correlation function (J(t w ), axis on the right). 

An other interesting quantity is I(t w ) = J d 3 xf(x,t w ) 
that is the equivalent of a susceptibility (it is equal to 
N({q(t w ,s t w ) 2 ) - (q(t w ,s t w )) 2 ). Roughly speaking it 
is proportional to the number of particles which move in 
a correlated way. If the typical event correspond to the 
formation of a fluidified domain of radius of order R(t w ), 
I(t w ) should be proportional to R(t w ) 3 . The correspond- 
ing exponent for / is about 0.52 (a similar conclusion was 
reached in [Q by an analysis restricted to much smaller 
samples, N < 258). It would be tempting to speculate 
that the exact exponents are 1/6 and 1/2. 

We have also similar data also for different values of 
s, i.e. s = 1.5 and s — .375. The dependence on s 
of R(t w ) is mild and the analysis of these data leads to 
similar conclusions. The data for f(r,t w ) at different s 
are indeed rather similar and there is mainly a difference 
in the absolute normalization. This fact is consistent 
with the possibility that the correlated movements of a 
large number of particles happen in a time interval which 
is much smaller of t w (this point should be studied more 
carefully). 

Our data indicate that the process of rearrangements 
of the systems happens on a spatial scale which becomes 
larger and larger as function of the time. Events that 
correspond to the rearrangements of regions of size R 
have a characteristic time which diverge approximately 
as R 6 . A comparison with the data of || shows that in 
our simulations we have explored a region of time where 
the relaxation of energy indicates that the approach to 
equilibrium is dominated by crossing of entropic barriers. 
For values of time near to the largest ones used in our 
simulations the energy relaxation enter in a new regime 
||, which is quite likely dominated by energetic barriers. 
It would be extremely interesting to extend these studies 
to longer times and larger systems, to see if there is any 
change in the behaviour of R(t w ) when we enter in this 



new regime. 

I thank Walter Kob and David Lancaster for useful 
discussions. 
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